GLOBAL ILL-POSEDNESS OF THE ISENTROPIC SYSTEM OF GAS 

DYNAMICS 

ELISABETTA CHIODAROLI, CAMILLO DE LELLIS AND ONDREJ KREML 

^y-v Abstract. We consider the isentropic compressible Euler system in 2 space dimensions 

T-H with pressure law p{p) = p^ and we show the existence of classical Riemann data, i.e. pure 

^^ jump discontinuities across a line, for which there are infinitely many admissible bounded 

^^ weak solutions (bounded away from the void) . We also show that some of these Riemann 

^~j data are generated by a 1-dimcnsional compression wave: our theorem leads therefore to 

, ^ Lipschitz initial data for which there are infinitely many global bounded admissible weak 

^^ solutions. 
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1. Introduction 



< 

r-| Consider the isentropic compressible Euler equations of gas dynamics in two space di- 

5^ mensions. This system consists of 3 scalar equations, which state the conservation of mass 

S and linear momentum. The unknowns are the density p and the velocity v. The resulting 

'— ' Cauchy problem takes the form: 

<^ r dtp + AiY^{pv) = 

en I dt{pv) + div^ {pv (^v) + V^[p(p)] = , . 

^ p(-,0) = p° ^^-^^ 

^ The pressure J9 is a function of p determined from the constitutive thermodynamic relations 

(-r^ of the gas under consideration and it is assumed to satisfy p' > (this hypothesis guarantees 

^ — I also the hyperbolicity of the system on the regions where p is positive). A common choice 

> is the polytropic pressure law p(p) = np^ with constants /t > and 7 > 1. The classical 

k> kinetic theory of gases predicts exponents 7 = 1 + |, where d is the degree of freedom 

^ of the molecule of the gas. Here we will be concerned mostly with the particular choice 

p{p) = p^. However several of our technical statements hold under the general assumption 
p' > and the specific choice p{p) = p^ is relevant only to some portions of our proofs. 

It is well-known that, even starting from extremely regular initial data, the system (1.1) 
develops singularities in finite time. In the mathematical literature a lot of effort has been 
devoted to understanding how solutions can be continued after the appearance of the first 
singularity, leading to a quite mature theory in one space dimension (we refer the reader 
to the monographs [1],[7] and [18]). In this paper we show that, in more than one space 
dimension, the most popular concept of an admissible solution fails to yield uniqueness even 
under very strong assumptions on the initial data. In particular we consider bounded weak 
solutions of (1.1), satisfying (1.1) in the usual distributional sense (we refer to Definition 
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3.1 for the precise formulation), and we call them admissible if they satisfy the following 
additional inequality in the sense of distibutions (usually called entropy inequality, although 
for the specific system (1.1) it is rather a weak form of energy balance): 

|2 



\v 
dt I pe{p) + p— - 1 + div^ 



P^ip) + P^r + pip) 



< (1.2) 



where the internal energy e : IR+ — ?> M is given through the law p{r) = r'^e'ir). Indeed, 
admissible solutions are required to satisfy a slightly stronger condition, i.e. a form of (1.2) 
which involves also the initial data, see Definition 3.2. For all solutions considered in this 
paper, p will always be bounded away from 0, i.e. p > Cq for some positive constant Cq. 
We denote the space variable as x = (xi,X2) € M^ and consider the special initial data 

( {p-,v^) if X2 <0 
{p%x)y{x)):=l (1.3) 

[ {p+,v+) iix2 > 0, 

where p±, v± are constants. It is well-known that for some special choices of these constants 
there are solutions of (1.1) which are rarefaction waves, i.e. self-similar solutions depending 
only on t and X2 which are locally Lipschitz for positive t and constant on lines emanating 
from the origin (see [7, Section 7.6] for the precise definition). Reversing their order (i.e. 
exchanging + and — ) the very same constants allow for a compression wave solution, i.e. 
a solution on M^x] — oo, 0[ which is locally Lipschitz and converges, for t ^ 0, to the jump 
discontinuity of (1.3). When this is the case we will then say that the data (1.3) are 
generated by a classical compression wave. 

We are now ready to state the main theorem of this paper 

Theorem 1.1. Assume p{p) = p^ . Then there are data as in (1.3) for which there are 
infinitely many hounded admissible solutions {p,v) of (1.1) on M? x [0, oo[ with infp > 0. 
Moreover, these data are generated by classical compression waves. 

It follows from the usual treatment of the 1-dimensional Riemann problem that for the 
data of Theorem 1.1 uniqueness holds if the admissible solutions are also required to be 
self-similar, i.e. of the form {p,v){x,t) = (r (^) ,w (^)) and to have locally bounded 
variation (see Proposition 8.1). Note that such solutions must be discontinuous, because 
the data of Theorem 1.1 are generated by compression waves. We in fact conjecture that 
this is the case for any initial data (1.3) allowing the nonuniqueness property of Theorem 
1.1: however this fact does not seem to follow from the usual weak-strong uniqueness (as 
for instance in [7, Theorem 5.3.1]) because the Lipschitz constant of the classical solution 
blows up as t 4 0- Related results in one space dimension are contained in the work of 
DiPerna [15] and in the works of Chen and Frid [3], [4]. 

As an obvious corollary of Theorem 1.1 we arrive at the following statement. 

Corollary 1.2. There are Lipschitz initial data (p°,f°) for which there are infinitely many 
bounded admissible solutions {p,v) of (1.1) on M^ x [0, oo[ with infp > 0. These solutions 
are all locally Lipschitz on a finite interval on which they all coincide with the unique 
classical solution. 
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We note in passing that, although the last statement of the corollary can be directly 
proved following the details of our construction, it is also a consequence of the admissibility 
condition, the Lipschitz regularity of the compression wave (before the singular time is 
reached) and the well-known weak-strong uniqueness of [7, Theorem 5.3.1]. 

1.1. /i-principle and the Euler equations. The proof of Theorem 1.1 relies heavily on 
the works of the second author and Laszlo Szekelyhidi, who in the paper [10] introduced 
methods from the theory of differential inclusions to explain the existence of compactly 
supported nontrivial weak solutions of the incompressible Euler equations (discovered in 
the pioneering work of Scheffer [19]; see also [20]). It was already observed by the same 
pair of authors that these methods could be applied to the compressible Euler equations 
and lead to the ill-posedness of bounded admissible solutions, see [11]. However, the data 
of [11] were extremely irregular and raised the question whether the ill-posedness was due 
to the irregularity of the data, rather than to the irregularity of the solution. 

A preliminary answer was provided in the work [5] where the first author showed that 
data with very regular densities but irregular velocities still allow for nonuniqueness of 
admissible solutions. The present paper gives a complete answer, since we show that even 
for some smooth initial data nonuniqueness of bounded admissible solutions arises after the 
first blow-up time. It remains however an open question how irregular such solutions have 
to be in order to display the pathological behaviour of Theorem 1.1. One could speculate 
that, in analogy to what has been shown recently for the incompressible Euler equations, 
even a "piecewise Holder regularity" might not be enough; see [13], [14], [16], [2] and in 
particular [8]. 

This paper draws also heavily from the work [22] where Szekelyhidi coupled the methods 
introduced in [10]-[11] with a clever construction to produce rather surprising irregular 
solutions of the incompressible Euler equations with vortex-sheet initial data. This work 
of Szekelyhidi was in turn motivated by the so-called Muskat problem (see [6], [23] and [21]; 
we moreover refer to [12] for a rather detailed survey). Indeed the basic idea of looking for 
piecewise constant subsolutions as defined in Section 3 stems out of several conversations 
with Szekelyhidi and have been inspired by a remark of Shnirelman upon the proof of [22]. 

1.2. Acknowledgements. The research of Camillo De Lellis has been supported by the 
SNF Grant 129812, whereas Ondfej Kreml's research has been financed by the SCIEX 
Project 11.152. The authors are also very thankful to Laszlo Szekelyhidi for several en- 
lightning conversations. 

2. Ideas of the proof and plan of the paper 

2.1. Subsolutions. Especially relevant for us is the appropriate notion of subsolution, 
which allows to use the methods of [10]- [11] to solve the equations and impose a certain 
specific initial data. We give here a brief description of the concept of subsolution relevant 
to us and refer to [12] for the motivation behind it and its links to existing literature in 
physics and mathematics. 
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Consider first some data as in (1.3). We then partition the upper half space {t > 0} in 
regions contained between half-planes meeting all at the line {t = X2 = 0} , see Definition 
3.3 and cf. Figure 1. We then define the density function p = p to be constant in each 
region: this density function will indeed give the final p for all the solutions we construct 
and it is therefore required to take the constant values p± in the outermost regions P±. 




Figure 1. A "fan partition" in five regions. 

We then solve the compressible Euler equations (1.1) in each region Pi, . . . ,P/v using 
the methods of [11]. Indeed observe that in each such region the density is constant 
and thus it suffices to to construct solutions of the incompressible Euler equations with 
constant pressure. Employing the methods of [11] we can also impose that the modulus 
of the velocity is constant (in each region): its square will be denoted by Cj. In [11] such 
solutions are constructed adding oscillations to an appropriate subsolution, which consists 
of a pair v,u of smooth functions, the first taking values in M^ and the second taking 
values in the space of symmetric, trace-free 2x2 matrices. These functions satisfy the 
linear system of PDEs 

dfV + div^M = 



divj^f = . 

and a suitable relaxation of the nonlinear constraints u = v ^ v — ■'^Id. 

In our particular case we will choose our subsolutions to be constant on each region 
Pf. the corresponding values will be denoted by {pi,Vi,Ui) and the corresponding glob- 
ally defined (piecewise constant) functions (p, v, u) will be called fan subsolutions of the 
compressible Euler equations. We then wish to choose our subsolution so that, after solv- 
ing (1.1) in each region Pi with the methods of [11], the resulting globally defined (p, f) 



GLOBAL ILL-POSEDNESS FOR COMPRESSIBLE EULER 5 

are admissible global solutions of (1.1). This leads to a suitable system of PDEs for the 
piecewise constant functions {p,v,u) which are summarized in the Definitions 3.4 and 3.5. 
In Section 3 we then briefly recall the notions of the papers [10]-[11] and in Section 4 we 
describe how to suitably modify the arguments there to reduce the proof of Theorem 1.1 
to the existence of the "fan subsolutions" of Definitions 3.4 and 3.5: the precise statement 
of this reduction is given in Proposition 3.6. 

2.2. The algebraic system. In Section 5, by making some specific choices, the existence 
of such subsolution is reduced to finding an array of real numbers satisfying some algebraic 
identities and inequalities, see Proposition 5.1. Indeed, since the functions (j),v,u) assume 
constant values in each region of the fan decomposition, these conditions are nothing but 
suitable "Rankine-Hugoniot type" identitites and inequalities. Although at this stage all 
computations can be carried in general, we restrict our attention to a fan decomposition 
which consists of only three regions. Therefore, the resulting solutions provided by Propo- 
sition 3.6 (and therefore also those of Theorem 1.1) will take the constant values {p±,v±) 
outside a "wedge" of the form Pi = {u^t < X2 < v+t}: inside this wedge the solutions will 
instead behave in a very chaotic way. 

Thus far, all the statements can be carried out for a general pressure law p. In the case 
p{p) = p^ we also compute explicitely the well-known conditions that must be imposed on 
the velocities v± and p± so that the corresponding data (1.3) are generated by a compression 
wave: this gives then an additional constraint. Observe that for such data the "classical 
solution" will be a simple shock wave traveling at a certain speed, whereas the nonstandard 
solutions of Theorem 1.1 "open up" the singularity and fill the corresponding region Pi 
with many oscillations. 

Coming back to the algebraic constrains of Proposition 5.1, although there seems to 
be a certain abundance of solutions to this set of identities and inequalities, currently we 
do not have an efficient and general method for finding them. We propose two possible 
ways in the Sections 6 and 7. That of Section 6 is the most effective and produces the 
initial data of Theorem 1.1 which are generated by a compression wave. That of Section 
7 is an alternative strategy, where, instead of making a precise choice of the pressure law 
p, we exploit it as an extra degree of freedom: as a result this method gives data as in 
Theorem 1.1 but with a different pressure law, which is essentially a suitable smoothing 
of the step-function. We also do not know whether any of these data are generated by 
compression waves. 

2.3. Classical Riemann problem. Finally in Section 8 we show that the self-similar 
solutions to (1.1)-(1.3) are unique: this follows from classical considerations but since we 
have not been able to find a precise reference, we include the argument for completeness. 

3. Subsolutions 

3.1. Weak and admissible solutions of (1.1). We recall here the usual definitions of 
weak and admissible solutions to (1.1). 
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Definition 3.1. By a weak solution of (1.1) on M^x [0, oo[ we mean a pair {p,v) G L°^(]R^x 
[0, oo[) such that the following identities hold for every test functions ip G C^(]R^ x [0, oo[), 



eC. 



00/10)2 



X [0,ooD: 



[pdt^p + pv ■ V xi^] dxdt + / p^{x)tlj{x,0)dx 



(3.1) 



'0 Jr^ 

/ / [pv ■ dt(l) + pv (g) v. Dx(t) + p{p) divx (/)] + p^{x)v^{x)-(j){x,0)dx = 0. (3.2) 

Jo Jr2 Jtr2 

Definition 3.2. A bounded weak solution {p,v) of (1.1) is admissible if it satisfies the 
following inequality for every nonnegative test function (/? G C^(]R^ x [0, oo[): 

I |2\ / I |2 

rU \ / \V\ 

P^ip) + P^ \dtV+ { P^ip) + P^ + P{p) 1 V ■ VxV 




p\x)eip'ix))+p%x) 



|t;0(a:)f 



ip{x, 0) dx > . 



(3.3) 



3.2. Subsolutions. To begin with, we state more precisely the definition of subsolution 
in our context. Here iSg^^ denotes the set of symmetric traceless 2x2 matrices and Id is 
the identity matrix. We first introduce a notion of good partition for the upper half-space 

M2x]0,oo[. 

Definition 3.3 (Fan partition). A fan partition of ]R^x]0, oo[ consists of finitely many 
open sets P_, Pi, . . . , P/v, P+ of the following form 

P_ = {{x,t) ■.t>0 and Xs < z/_t} (3.4) 

P+ = {(a;,t) : t > and X2 > u+t} (3.5) 

Pi = {{x,t) : t > and I'i-it < X2 < z/jt} (3.6) 

where z/_ = z/q < z^i < . . . < z/tv = z^+ is an arbitrary collection of real numbers. 

The next two definitions are then motivated by the discussion of Section 2.1. However at 
the present stage it is not completely clear why the relevant partial differential equations 
(and inequalities!) for the piecewise constant solutions are given by (3.8), (3.9) and (3.10): 
their role will become transparent in the next subsection when we prove Proposition 3.6. 

Definition 3.4 {Fan Compressible subsolutions). A fan subsolution to the compressible 
Euler equations (1.1) with initial data (1.3) is a triple {p,v,u) : ]R^x]0,oo[-> (]R+,]R2,>So''^) 
of piecewise constant functions satisfying the following requirements. 

(i) There is a fan partition P_, Pi, . . . , P/v, P+ of IR^xjO, oo[ such that 

N 
{p,V,u) = ^{pi,Vi,Ui)lp. + {p_,V_,U_)lp_ + {p+,V+,U+)lp^ 

1=1 
where pi, Vi, Ui are constants with pi > and u± = v±<g v± — ^|f±pld; 
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(ii) For every i G {1, . . . , N} there exists a positive constant Cj such that 

C- 
Vi®Vi-Ui < yid. (3.7) 

(iii) The triple (p, v, u) solves the following system in the sense of distributions: 

9ip + div,(pt;) = (3.8) 

dt{-pv) + div^ {-pu) + V^ I p(p) + 2 Zl CiP^'^P^ + Pl^^Plp+uP- I I = (3-9) 

Definition 3.5 (Admissible fan subsolutions) . A fan subsolution {p,v,u) is said to be 
admissible if it satisfies the following inequality in the sense of distributions 

Lp+uP- 
.... .... v^ / V^ 

dt (pi y IpJ + div. ( P. ^ ^ IpJ < . (3.10) 



dt ipe{p)) + div^ [ipe{p) + p{p)) v] + dt [ p^— lp+up_ + div^ P^^vl 



+5: 



2 ^^+^^- 


+ divj: 
/ 


V 2 




< 0. 





It is possible to generalize these notions in several directions, e.g. allowing partitions with 
more general open sets and functions Vi,Ui and pi which vary (for instance continuously) 
in each element of the partition. It is not difficult to extend the conclusions of the next 
subsection to such settings. However we have chosen to keep the definitions to the minimum 
needed for our proof of Theorem 1.1. 

3.3. Reduction to admissible fan subsolutions. Using the techniques introduced in 
[10]- [11] we then reduce Theorem 1.1 to finding an admissible fan subsolution. The precise 
statement is given in the following proposition. 

Proposition 3.6. Let p be any C^ function and {p±,v±) be such that there exists at least 
one admissible fan subsolution {p,v,u) of (1.1) with initial data (1.3). Then there are 
infinitely many bounded admissible solutions (p, f) to (1.1)-(1.3) such that p = p. 

The core of the proof is in fact a corresponding statement for subsolutions of the incom- 
pressible Euler equations which is essentially contained in the proofs of [10]- [11]. However, 
since our assumptions and conlusions are slightly different, we state them in the next 
lemma. 

Lemma 3.7. Let {v,u) G M^ x S^^^ and C > be such that v ®v — u < |-Id. For any 
open set f] C M^ x M there are infinitely many maps {v_,u) G L°^(]R^ x M, M^ x iSg^^) with 
the following property 

(i) y_ and u vanish identically outside Vt; 
(ii) diVa; V = and dtv + div^ u = 0; 
(iii) (v + v) {v + v) — {u + u) = |-Id a. e. on ^2. 
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The proof is a minor variant of the ones given in [10]- [1 1] but since none of the statements 
present in the hterature matches exactly the one of Lemma 3.7 we give some of the details 
in the next Section, referring to precise lemmas in the papers [10]-[11]. For the moment 
we show how Proposition 3.6 derives from Lemma 3.7. 



Proof of Proposition 3.6. We apply Lemma 3.7 in each region Vt 
any pair of maps given by such Lemma. Hence we set 



Pi and we call {Vi,Ui} 



N 



i=l 

N 



V := V + y ^Vj 

i=\ 

N 

u := u+ }^Ui 

i=l 



(3.11) 
(3.12) 



whereas p = p (as claimed in the statement of the Proposition!). We next show that the 
pair {p,v) is an admissible weak solution of (1.1)-(1.3). First observe that diVx{piVj) = 
since pi is a constant. But since v^ is supported in Pi and p = p = pi on Pi, we then 
conclude div x{pv.i) = 0. Thus we have 

dtp + div.x{pv) = dtp + div^ ipv + y^^pVi] 
=dtp + divx{pv) + ^ div x{pVi) = dtp + div^(pl;) = (3.13) 

i 

in the sense of distributions. Moreover, observe that 



V (^ V 



on P4 
on P_ 

{vi + Vi) (g) {vi + Vi) =Ui + Ui + ^Id on Pi 



V- (g V- 



and 



t;+ «+ — l|f+pld on P+ 



u= < 



f _ (X) V- 



Ui 



l|f_Pld on P_ 



on P . 



Moreover, on each region Pj we have 



Cip, 



\v\\\ pC,^^ _ 
pv(^v = p\v^v —id I + ^^Id = pu + piUi + 



Id. 
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Hence, we can write 



dt{pv) + diY^{pv (^v) + Vx[p{p)] = dtipv + ^ piV^ J + div^ I P^ + X^ PiVa ) 

--dt (pv) + div^ (pu) + Va; I p{p) + 2 X] CiPdp, + ^ 
+ y^ Pi dtVi + div^ Ui . 



gl^-lVlp- + 2l^+lV+lp+ 



f_|V-lP- + ol^+lV+lp+ 



(3.14) 



Therefore, by Definition 3.4 we conclude dt{pv) + div^(pf ®v) + V x[v{.p)] = 0. 
Next, we compute 

dt I pe{p) + ^p j + div^ ( ( pe{p) + ^p + p{p) j v 



--dt f pe{p) + Y^ -CiPilp^ + 



12 I 12 

-^P-lp_ + ^-P+Ip+ 



+ div^ 



pe{p) + p{p) + Y^ -CiPilp^ + 



— p_lp_ + ^-p+lp^ ) I V 



Using the condition (3.10) we therefore conclude 

dt [pe{p) + —p\ + div^ ( [pe{p) + — p + p{p) ] v 



1 



< ^div^ vApe{p)+p{p) + Y T^Cipilp, + — — p_lp_ + ^— p+lp+ 



+ E^^ 



(3.15) 



in the sense of distributions. Observe however that the function g is constant on each Pj, 
on which v^ is supported. Thus 



dt ( P^(p) + ^^p) +div^ ( (p£(p) + ^P + P(p)l ^1 < J^^div 



xVi=0. 



(3.16) 



So far we have shown that (3.1), (3.2) and (3.3) hold whenever the corresponding test 
functions are supported in M^xjO, cxd[. However observe that, since as r J, the Lebesgue 
measure of Pj fl {t = r} converges to 0, the maps p(-, r) and v{-, r) converge to the maps 
p° and f° of (1.3) strongly in Lj^^. This easily implies (3.1), (3.2) and (3.3) in their full 
generality. For instance, assume ip G C^(]R^x] — oo,oo[) and consider a smooth cut-off 
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function i} of time only which vanishes identically on ] — oo, e] and equals 1 on ]5, oo[, where 
< e < 6. We know therefore that (3.1) holds for the test function ipi), which implies that 

■& [pdtip + pv ■ Vxi^] dxdt + / / ■d'{t)p{x, t)il){x, t)dxdt = . 
'o Jw? Jo Jr'2 

Fix S and choose a sequence of i) converging uniformly to the function 

f if t < 

r]{t) = I 1 iit>6 

[ I ii0<t<6 

and such that their derivatives -i?' converge pointwise to |l]o,<5[- We then conclude 



T] [pdtip + pv ■ Vx^] dxdt + T / / pix, t)^lj{x, t)dxdt = . 

JR^ (^ Jo iR2 

Letting 6 10 we conclude (3.1). 

The remaining conditions (3.2) and (3.3) are achieved with analogous arguments, which 
we leave to the reader. D 

4. Proof of Lemma 3.7 

4.1. Functional set-up. We define Xq to be the space of {v, u) G C^{il, M^ x iSq^^) which 
satisfy (ii) and the pointwise inequality {v + v)<^{v + v) — {u + u) < |-Id. We then take the 
closure X of Xq in the L°° weak* topology and recall that, since X is a bounded (weakly*) 
closed subset of L°^ such topology is metrizable on X, giving a complete metric space 
{X, d). Observe that any element in X satisfies (i) and (ii) and we want to show that on a 
residual set (in the sense of Baire category) (iii) holds. We then define for any N eN\ {0} 
the map In as follows: to {v, u) we associate the corresponding restrictions of these maps 
to Bn{0)x] — N,N[. We then consider I^ as a map from (X, d) to F, where Y is the space 
L°°{Bn{0)x] - N, A^[,]R^ X ^0^^) endowed with the strong L"^ topology. Arguing as in [10, 
Lemma 4.5] it is easily seen that I^ is a Baire-1 map and hence, from a classical theorem 
in Baire category, its points of continuity are a residual set in X. We claim that 
(Con) if {v,u) is a point of continuity of In, then (iii) holds a.e. on Bn{0)x] — N, N[. 
(Con) implies then (iii) for those maps at which all In are continuous (which is also a 
residual set). 

The proof of (Con) is achieved as in [10, Lemma 4.6] showing that: 

(CI) If {v,u) G Xq, then there is a sequence {vk,Uk) C Xq converging weakly* to {v,u) 
for which 

2 



liminf \\v + Vk\\L2(r) > \\v + ^||L2(r) + /^ [C\^\ 



\v + v\\l2(^r-) 



where F = Bn{0)x] — N, N[ and /3 depends only on F. 

Indeed assuming that (CI) holds, fix then a point {v, u) E X where In is continuous 
and assume by contradiction that (iii) does not hold on F. By definition of X there is a 
sequence (w^, Uf.) C Xq converging weakly* to {v, u). Since the latter is a point of continuity 
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for Jat, we then have that v,.^v strongly in L'^iT). We apply (CI) to each {vi.,Ui.) and 
find a sequence {{vkj,Ukj)} such that 

2 



liminf \\v + VkjWL^r) > \\v + Vk\\h{T) + P {c\^\ - \\v 



^A:llL2{r) 



and (ffcjjMfcj) ^* {v.ki1kk)- ^ standard diagonal argument then allows to conclude the 
existence of a sequence {vk,j(k),Ukj(k)) which converges weakly* to {v,u) and such that 

limmi \\v + Vk,j(k)\\LHr) > Il^' + ^lli2(r) + P [C\r\ - ||t^ + t^||^2(r)j > ||^ + ^!||i2(r) • 

However this contradicts the assumption that {v,u) is a point of continuity for In- 

In order to construct the sequence of (CI) we appeal to the following Proposition and 
Lemma. 

Proposition 4.1 (Localized plane waves). Consider a segment a = [—p,p] C M^ x iSg^^, 
where p = A[(a, a §>> a) — (&, 6 (X> b)] for some A > and a ^ ±b with \a\ = \b\ = \fC . Then 
there exists a pair {v,u) G C^(-Bi(0)x] — 1, 1[) which solves 

dfV + diVxU = 

(4.1) 
divj-f = 

and such that 

(i) The image of {v,u) is contained in an e -neighborhood of a and J{v,u) dxdt = 0; 
(ii) J \v{x,t)\ dxdt > aX\b — a\, where a is a positive constant depending only on C. 

In order to the state the next lemma, it is convenient to introduce the following notation. 

Definition 4.2. Let C > be the positive constant of Lemma 3.7. We let U be the subset 
of M^ X (Sq^^ consisting of those pairs (a, A) such that a® a — A < ^Id. 

Lemma 4.3 (Geometric lemma). There exists a geometric constant Cq with the following 
property. Assume {a., A) G lA. Then there is a segment a as in Proposition ^.i with 
{a,A)-\-aCU and X\b — a\ > Co(C — lap). 

We are now ready to prove (CI). Let {v,u) G Xq. Consider any point {xQ,to) G F and 
observe that {v,u) + {v,u) takes values in U. Let therefore a be as in Lemma 4.3 when 
{a, A) = {v,u) + {v{xo,to),u{xo,to)) and choose r > so that {v,u) + {v{x,t),u{x,t)) + 
0" C W for any {x,t) G Br{xo)x]tQ — r, to + r[. For any e > consider a pair {v,u) 
as in Proposition 4.1 and define {Vxf^^to,r,Uxo,to,r){x,t) := {v,u) (^^^, ^-jr-)- Observe that 
{v.,u) + ('^xo,to,r?'"xo,io,r) £ -^0 provided e is sufficiently small, and moreover 

■Jxo,to,r\ >CoaX{C-\v + v{xo,to)\^y. (4.2) 

By continuity there exists Tq such that the conclusion above holds for every r < Tq and 
every (x, t) with Br{x)x]t — r, t + r[c F. Fix now k E N with | < rg. Set r := ^ and 
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find a finite number of points {xj,tj) such that the sets Br{xj)x]tj — r,tj + r[ are pairwise 
disjoint, contained in F and satisfy 

J2{C-\v + v{xj,tj)\^)r^>c(c\T\- I \v + v{x,t)\^dxdt\ , (4.3) 

where c is a suitable geometric constant. We then define 

J 

Since the supports of the {vx-^t,r,Ux.^f,r) are pairwise disjoint, {vk,Uk) belongs to Xq as 
well. Moreover, using the property that f{vx^t,r,Ux-,t,r) = 0, it is immediate to check 
that {vk,Uk) ^* {v.,u) in L'^- On the other hand it also follows from (4.2) and (4.3) that 



\Vk -miL^ 



ir)>ci(c\T\- f\v + vA 



where the constant Ci is only geometric. Using the weak* convergence of {vk,Uk) to {v,u) 
we can then conclude 



liminf lit; + t;fc 11^2 /r) = 11^ + v\\l2(Tp) + liminf ||t;fc - vW"^ 

k ' ' k 

> \\v + ^|li2(r) + |r| fliminf ||wfc - v\\lA 



> \\v + vWhiD + clW^W (c|r|- / \v + v\^^ 

which concludes the proof of the claim (CI). 

4.2. Proof of Proposition 4.1 and of Lemma 4.3. 

Proof of Proposition 4-1- Consider the 3x3 matrices 

Apply [11, Proposition 4] with n = 2 to Ua and Ub and let A{d) be the corresponding 
linear differential operator and r^ G M^ x Mt the corresponding vector. Let (p he a. cut-off 
function which is identically equal to 1 in -81/2(0) x] — ^,^[, is compactly supported in 
Bi{0)x] — 1, 1[ and takes values in [—1, 1]. For N very large, whose choice will be specified 
later, we consider the function 

0(x, t) = —\N~^ sm(Nr] ■ (x, t))ip{x, t) =: k{x, t)ip{x, t) 

and we let U{x,t) := A(9)(0). According to |11, Proposition 4], f/ : M^ x M — )■ S^^^ is 
divergence free and trace-free and moreover f/33 = 0. Note also that f^ (o)xi-i if ^(^) ^) — ^■ 
Define 
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It then follows easily that {v,u) satisfies (4.1) and that it is supported in Bi{0)x] — 1, 1[. 
Also, since A{d) is a 3rd order homogeneous linear differential operator with constant 
coefficients, \\U — ipA{d){K,)\\o < CXN~^, where C depends only on the cut-off function (/?: 
in particular we can assume that ||f/ — (/9A(9)(fi:)||o < e. On the other hand [11, Proposition 
4] clearly implies that 

we therefore conclude that U takes values in an e-neighborhood of the segment [—X{Ua — 
Ub),\{Ua — Ub)]- This obviously implies that {v,u) takes values in an e-neighborhood of 
the segment a. Finally, Let By2 be the 3-dimensional space-time ball in M^ x M, centered 
at and with radius |. Observe that 



■°l/2 '^-°l/2 



\v{x,t)\> I X\a — b\\cos{N{x,t) ■ ri)\dxdt = X\a — b\ j \cos{Nt\ri\)\dxdt . 
Moreover, 



lim / \cos{Nt\'r]\)\dxdt = a 



TVtoo / o3 

for some positive geometric constant a. D 

Proof of Lemma 4 ■ 3. Consider the set 

K^ := I {v, m) e M^ X S^""^ ■.u = v®v--ld, |i;p = c\ ■ 

It then follows from [11, Lemma 3] that U is the interior of the convex hull of K^. The 
existence of the claimed segment a is then a corollary of [11, Lemma 6], since A|6 — a| is 
indeed comparable (up to a geometric constant) to the length of a. D 

5. A SET OF ALGEBRAIC IDENTITIES AND INEQUALITIES 

In this paper we actually look at fan subsolutions with a fan partition consisting of only 
three sets, namely P-,Pi and P+. 

We introduce therefore the real numbers a, f3,'y,S,V-i,V-2,v+i,v+2 such hat 

(5.1) 
(5.2) 
(5.3) 

(5.4) 

Proposition 5.1. Let N = 1 and P_,Pi,P+ be a fan partition as in Definition 3.3. The 
constants fi,f_,f+,-Ui,p_,p_|_,pi as in (5.1)-(5.4) define an admissible fan subsolution as 
in Definitions 3.4-3.5 if and only if the following identities and inequalities hold: 



Vl 


= («,/?), 


V- 


= {v^i,V^2) 


v+ 


= {v+l,V+2) 


Ml 


_ / 7 6 
~ I 5 -7 
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• Rankine-Hugoniot conditions on the left interface: 

u_{p_-pi) = p_v_2-pi(3 (5.5) 

P-{p-V-i - pia) = p-V-iV-2 - pi6 (5.6) 

C 

i/_(p_t;_2-pi/3) = p-t;!2 + Pi7 + p(P-) -p(Pi) -Piy; (5.7) 

• Rankine-Hugoniot conditions on the right interface: 

^+{pi - P+) = Pi/3 - P+v+2 (5.8) 

z/+(pia - p+f+i) = pi5 - p+Vj^iv+2 (5.9) 

C 

v+{piP - P+v+2) = -Pi7 - P+^+2 + P(Pi) - P(P+) + P^Y' ^^'^'^^ 

• Suhsolution condition: 

a^ + P'^<Ci (5.11) 



2 -« +7) ( Y-/3 -7)-(5-«/3)'>0; (5.12) 

Admissibility condition on the left interface: 

I P r" 

Z/_(p_£(p_) - pi£(pi)) + V- I p_^ Ply 

< [(p-e(p-) + p(p-))t;-2 - (Pi^(Pi) + p(Pi))/3] + I P-^-2^ - Pi/3^ 1 ; (5.13) 



2 ^ ' 2 



Admissibility condition on the right interface: 



Ci \v. 



|2 



«^+(Pl^(Pl) - P+^{P+)) + '^+\ Pl^ - P+ 



2 '^^ 2 
< [(Pi^(Pi) + P(Pi))/3 - iP+<P+) + P(P+))^+2] + I Pi/3 Y - P+vJ-^ I . (5.14) 



Proof. Observe that the triple (p, v, u) does not depend on the variable xi. We will therefore 
consider it as a map defined on the t, X2 plane. The various conditions and inequalities 
follow from straightforward computations, recalling that the maps p, v and u are constant 
in the regions P_, Pi and P+ shown in Figure 2. In particular 

• The identities (5.5) and (5.8) are equivalent to the continuity equation (3.8), in 
particular they derive from the corresponding "Rankine-Hugoniot" type conditions 
at the interfaces between P~ and Pi (the left interface) and Pi and P+ (the right 
interface), respectively. 
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Figure 2. The fan partition in three regions. 

The identitities (5.6) and (5.9) are the Rankine-Hugoniot conditions at the left and 

right interfaces resulting from the first component of the momentum equation (3.9); 

similarly (5.7) and (5.10) correspond to the Rankine-Hugoniot conditions at the left 

and right interfaces for the second component of the momentum equation (3.9). 

The inequalities (5.11) and (5.12) are derived applying the usual criterion that the 

matrix 

C 

M := —\d-vi®vi+ui (5.15) 

is positive definite if and only if tr M and det M are both positive. 
Finally, the conditions (5.13) and (5.14) derive from the admissibility condition 
(3.10), again considering, respectively, the corresponding inequalities at the left 
and right interfaces. 

D 



6. First method: data generated by compression waves for p{p) = p^ 

In this section we show how to find solutions of the algebraic constraints in Proposition 
5.1 when p(p) = p^ with pairs (p±,f±) which can be connected by a compression wave, 
thereby showing Theorem 1.1. We start by recalling the following fact, which can be 
easily derived using (by now) standard theory of hyperbolic conservation laws in one space 
dimension. 



Lemma 6.1. Let < p < p^ 



( — ^,0) andv^ = (--^. 



.2V2{,/p^- TpT))- Then 
there is a pair {p,v) E W^^^ D L°°(M?x] - oo,0[,]R+ x M^) gu^h that 



(i) P+>P>P-> 0; 
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(ii) The pair solves the hyperbolic system 

( dtp + div^ipv) = . . 

with p{p) = p^ in the classical sense (pointwise a.e. and distrihutionally) ; 
(iii) /or 1 1 the pair {p{-,t),v {■,{)) converges pointwise a.e. to (p°,f°) as in (1.3); 
(iv) (p(-, t),v{-, t)) e W^^°° for every t < 0. 

As already mentioned, the proof is a very standard application of the one-dimensional 
theory for the so-called Riemann problem. However, we give the details for the reader's 
convenience. 

Proof. We look for solutions (p, v) with the claimed properties which are independent of 
the Xi variable. Moreover we observe that, since we will produce classical Wf^^ solutions, 
the admissibility condition (3.3) will be automatically satisfied as an equality because 

I 12 / I 12 

\v\ I , . \v\ 



p^{p) + ^p> {p^{p) + -yp + P^P) ) ^ 

is an entropy-entropy flux pair for the system (6.1) (cf. [7, Sections 3.2, 3.3.6, 4.1]). We 
then introduce the unknowns 

{mi{x2,t),m2{x2,t)) = m{x2,t) := v{x2,t)p{x2,t) 
and hence rewrite the system as 

dtp + dx2^2 = 

dtm^ + d,, hf^) = (6.2) 

9im2 + 9,, (^ + p2) = 

Observe that if {p,m) is a solution of (6.2) then so is 

(p(x2,t),m(a;2,t)) := {p{-X2,-t),m{-X2,-t)) . 

Moreover, if (p, m) is locally Lipschitz and hence satisfies the admissibility condition with 
equality, so does {p,rh). We have therefore reduced ourselves to finding classical Wi 
solutions on ]Rx]0, oo[ of (6.2) with initial data 

if X2 < 0, 



l,oo 
loc 



Po(:r) :=<;-- ^^^^^^JJ' (6.3) 



and 

^^^^^_^^Un:^{-f^,2V-2pniVp-L-Vp-n)) ^f -^ > 0, ^^^^ 

[mL:= (-1,0) iix2<0, 

where p^ = pi > Pr = p~ > 0, rriL = v^p+ and rriR = f_p_. 

The problem amounts now in showing that, under our assumptions, there is a classical 
rarefaction wave solving, forward in time, the system (6.2) with initial data (po,'n2o) as in 
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(6.3) and (6.4). We set therefore p{p) = p^ and we look for a locally Lipschitz self-similar 
solution {PjItl) to the Riemann problem (6.2)-(6.3)-(6.4): 

{p,m){x2,t) = {R,M) ( — j , -oo < a;2 < cx), < t < oo . (6.5) 

Thus (-R, M) are locally Lipschitz functions on (— oo, oo) which satisfy the ordinary differ- 
ential equations 

d 



^^[M,{O-mO]+R{O = 



d 

d 



^^^fP-eM.(0' 



+ Ml (0=0 



+ p(^(0)-eM2(0 



M2(O = 0. 



L R{0 

Before analyzing our specific Riemann problem, we review some general notions for system 
(6.2) (referring the reader to the monographs [7] and [18]). If we introduce the state vector 
U := (p, 7721,7722), we can recast the system (6.2) in the general form 

dtU + d^,F{U) = 0, 

where 




FiU) :- 
By definition (cf. [7]) the system (6.2) is hyperbolic since the Jacobian matrix DF{U) 

DF{U) = 





-mini2 



+ P'ip) 




has real eigenvalues 



X. = '-^-v9{p), A2 = ^, M='-^ + V¥{p) (6.6) 

P P P 



and 3 linearly independent eigenvectors 

1 




R,= \ 7 \, R2=\ I \, Rs=\ 7 I • (6.7) 



7 - vp'ip) J V / V ? + vp'(p) 

The eigenvalue Aj of DF^ i = 1,2, 3, is called the i- characteristic speed of the system (6.2). 
On the part of the state space of our interest, with p > 0, the system (6.2) is indeed strictly 
hyperbolic. Finally, one can easily verify that the functions 



7772 



p 



V(^) , _ 777i _ 7772 f \/¥ij) 



W3 = h / — dr, W2 = — , wi = / — dr (Q.i 

P Jo ^ P P Jo T 
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are, respectively, (1- and 2-), (1- and 3-), (2- and 3-) Riemann invariants of the system 
(6.2) (for the relevant definitions see [7]). 

In order to characterize rarefaction waves of the reduced system (6.2), we can refer to 
Theorem 7.6.6 from [7]: every i- Riemann invariant is constant along any i- rarefaction 
wave curve of the system (6.2) and conversely the i- rarefaction wave curve, through a 
state (p, m) of genuine nonlinearity of the i- characteristic family, is determined implicitly 
by the system of equations Wi{p,m) = Wi(j),fn) for every i- Riemann invariant Wi. As an 
application of this theorem, we obtain that (p^, m/j) lies on the 1- rarefaction wave through 
(Plj'^^l)- Indeed, the 1- rarefaction wave of the system (6.2) through the point (pi,mi) 
is determined in terms of the Riemann invariants ^3 and W2 by the equations 



mi = , m2 = p ttT, (6.9) 

PL Jp T 

with p < pl. In the case of pressure law p{p) = p^, the equations (6.9) read as 

mi = -^, m2 = 2v^p (V^ - ^) . (6.10) 

Pl 

Clearly, the constant state {pu, rriR), as defined by (6.3)-(6.4), satisfies the equations (6.10). 

Since, according to Theorem 7.6.5 in [7], there exists a unique 1-rarefaction wave through 

(Plj'/til), we have shown the existence of our desired self-similar locally Lipschitz solution. 

Observe that, by construction, p+ = Pl > p > p_r = p_ > 0, thereby showing (i). The 

claim (iv) follows easily because there exists a constant C > such that, for every positive 

time t, the pair {p.,m) takes the constant value {pR,Tnfi) for X2 > Ct and {pL,mL) for 

X2 < -Ct. D 

We next show the existence of a solution of the algebraic constraints of Proposition 5.1 
such that in addition (p±,f±) satisfy the identities of Lemma 6.1. 

Lemma 6.2. Let p{p) = p^ . There exist p±,v± satisfying the assumptions of Lemma 6.1 
and pi,Ci,vi,ui,i'± satisfying the algebraic identities and inequalities (5.5)-(5.14). 

Proof. Taking into account that p(p) = p^ and therefore e{p) = p, we substitute the 
identities of Lemma 6.1 into the unknowns of Proposition 5.1 and reduce (5.8)-(5.10) to 

^+{pi - P+) = PiP (6-11) 

u+{pia + 1) = pi5 (6.12) 

C 
Z/+pi/3 = -P17 + P1 -p^ + piy. (6.13) 

Similarly, we reduce (5.5)-(5.7) to 

^-(P- - Pl) = 2V2p_(v/pT - ypT) - pi/3 (6.14) 

^- " "'- -2y2^(v/pT-v^)-pi5 (6.15) 



iy_[ pia) 

P+ P^ 

z/_(2v^p_(yp7 - ^) - pi/3) = 8p_(v/p7 - v^)' + Pi7 + P- - P? - Pi^ • (6.16) 



2 
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The identities of Lemma 6.1 do not influence the form of (5.11)-(5.12). Instead, plugging 
them into (5.13)-(5.14) the latter are reduced to 



z/_ \p_ 

< V2p-{ 



P\ 



P 



1 



/p-; 



Cipi 



P+ 



'P-> 



1p\fi- 



(^Cipi 



,2 2 I CiPi 



1 



2PH 



<2pl(3 



Cip.P 



(6.17) 
(6.18) 



We next make the choice z/+ = /3 = 5 = and hence (6.11), (6.12) and (6.18) are automat- 
ically satisfied. The remaining constraints above then become 







Pi7 + P?-P+ + Pi- 



z/_(p_-pi)=2y2p_ 
P- 



^Pa 



'P-] 



Pittj 



2^2^ 
P^ 



z/_ 



(2v^P_(Vp7 - VpT)) =8p_(v/pT - Vp^f + Pi7 + P- - P? - Pi 



Ci 



(6.19) 
(6.20) 
(6.21) 

(6.22) 



and 



z/_ p 



P\ 



P- 



2p\ 



4p- 



<v^p_(v/pT-v^)(4p- 



/P+ 
1 



Cipi 



+ 8p-(yp^ 



^p-; 



Moreover, (5.11) and (5.12) become 

c? <Ci 



<(^-a2 + 7 



2 



7 



(6.23) 

(6.24) 
(6.25) 



Summarizing we are looking for real numbers v_ < 0, < p_ < p+, pi, a, 7 and C\ satisfy- 
ing the set of identities and inequalities (6.19)-(6.25). 

We next choose p_ = 1 < 4 = p_(_ and simplify further (6.19)-(6.23) as 

CiPi , „2 „., .. ^0 ^Q26) 

(6.27) 
(6.28) 

2v^z/_ (6.29) 



2 


-^Pi- 


P\l - i"^ 


^-(: 


L-Pi) = 


2v^ 


z/_ 


- + api 


/ " 2 


9 + 


Pi7 - p\ 


Cipi 
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-(-^-?-^)^^(-^^) 



(6.30) 



T_ 

^ Uii^ y^j.^Kj, yKj.^^j iiiiJJiJ ^- 2^- 

Therefore, our constraints further simphfy to looking for pi,7, Ci such that 



We now observe that (6.27) and (6.28) imply a = --\ and (6.26)-(6.29) imply z/. 

;h 
1 



^Q < C, (6.31) 

0<(^-l..)(^-.) (6.32) 

= ^ + P?-pi7-16 (6.33) 

8 = -7(l-pi) (6.34) 

48 + i>-7(5.1-,;-%^)^ (6,35) 

From (6.34) we derive pi = y and inserting this into (6.33) we infer ^ — 7 = f§|- In turn 
this last identity reduces (6.32) to the inequality 

C, > 1 + ^ . (6.36) 

16 105 ^ ' 

The remaining constraints (6.33) and (6.35) simplify to: 

Y-^=i05 ('-'^^ 

48+i + 35 + |-f >lfl. (6.38) 

We therefore see that 7 can be obtained from Ci through (6.37). Hence the existence of 
the desired solution is equivalent to the existence of a Ci satisfying (6.36) and (6.38). Such 
Ci exists if and only if 

15 / 1 559\ ,1 7 225 

y(l6 + I05J*'"^+i+^^+32--- 
which can be trivially checked to hold. D 

Theorem 1.1 and Corollary 1.2 now easily follow. 

Proofs of Theorem 1.1 and Corollary 1.2. Let p{p) = p^ and consider the p±.,v± given by 
Lemma 6.2. Applying Propositions 5.1 and 3.6 we know that there are infinitely many 
admissible solutions of (1.1)- (1.3) as claimed in the Theorem. 

Let now {pf,Vf) be any such solution and let {pb,Vb) be the locally Lipschitz solutions 
of (1.1) given by Lemma 6.1. It is straightforward to check that, if we define 

{Pf,vf){x,t) ift>0 

{p,v){x,t) := { (6.39) 

{pb,Vb){x,t) ift<0. 
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then the pair (p,f) is a bounded admissible solution of (1.1) on the entire space-time 
M^ X M with density bounded away from 0. Moreover (p(-, t),v{-, t)) is a bounded Lipschitz 
function for every t < 0. In particular we can define (p, v){x, t) = (p, v){x, t — 1) and observe 
that, no matter which of the infinitely many solutions {pf,Vf) given by Theorem 1.1 we 
choose, the corresponding {p,v) defined above is an admissible solution as in Corollary 1.2 
for the bounded and Lipschitz initial data {p^,v^) = {pi),Vb){-, —1). □ 

Remark 6.3. In fact, it is not difficult to see by a simple continuity argument that the 
conclusions of Theorem 1.1 and Corollary 1.2 hold even for general pressure laws p{p) = p"' 
with 7 in some neighborhood of 2. 

7. Second method: further Riemann data for different pressures 

In this section we describe a second method for producing solutions to the algebraic set 
of equations and inequalities of Proposition 5.1. Unlike the method given in the previous 
section, we do not know whether this one produces Riemann data generated by a compres- 
sion wave: we can only show that this is not the case for the ones which we have computed 
explicitely. Moreover we do not fix the pressure law but we exploit it as an extra degree 
of freedom. On the other hand the reader can easily check that the method below gives a 
rather large set of solutions (i.e. open) compared to the one of Lemma 6.2 (where we do 
not know whether one can perturb the choice z/+ = 0). 

Lemma 7.1. Set v± = (±1,0). Then there exist v±, p±, pi, a,P,'-f,S, Ci and a smooth 
pressure p with p' > for which the algebraic identities and inequalities (5.5) -(5. 14) are 
satisfied. 

7.1. Part I of the proof of Lemma 7.1: reduction of the admissibility conditions. 

We rewrite the conditions (5.5)-(5.10) 

z/_(pi-p_) =pi/3 (7.1) 

u-ip- + pia) = pi5 (7.2) 

C 

Piy -Pi7 + p(Pi) -P(P-) = t^^Pi(3 (7.3) 

Mpi - P+) = Pi^ (7.4) 

z/+(pia-p+) = pi6 (7.5) 

C 
Pi^ - Pil + P{pi) - P{p+) = v+piP- (7.6) 

The conditions (5.11) and (5.12) are not affected by our choice. The conditions (5.13) and 
(5.14) become 

^- ip-^p-) - Pi^pi) + Y - /^i y) + /^ (pi^Pi) + P{pi) + Pi y) ^0 (7-7) 

i^+ (pi^(Pi) - P+e{p+) + Pi Y - y) " '^ (pi^(Pi) + P(Pi) + Pi y) - ° • ('^•^^ 
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Plugging (7.1) and (7.4) into, respectively, (7.7) and (7.8) we obtain 

^+ \pie{pi) - P+e{p+) + Pi^^^^Y) ~ ^ y^^^P'^ + P^P^^ + '^^^^^) - ° • ^'^•^°^ 

We next rely on the following 

Lemma 7.2. Assume that 

z/_<0<z/+, (7.11) 

P-<P+- (7.12) 

Then, there exist pressure functions p G C°^([0, +oo[) with p' > on]0, +oo[ such that the 
admissibility conditions (7.9) -(7. 10) for a subsolution are implied by the following system 
of inequalities: 

C — 1 
(P(P+) - Pipi)) (P+ - Pi) > ^ — P+Pi (7.13) 

(p(Pi) - Pip^)) (Pi -p-)> ^^P^Pi- (7.14) 

Proof. First, let us define g{p) := pe{p). In view of the relation p(p) = p'^e'{p), we obtain 

^ p = ^ p + — • 
p 

Thus, by vitue of (7.1) and (7.4), respectively, we can rewrite (7.9) and (7.10) as follows: 

u-{g{p-) - g{p^)) + z/_(Pi - P-)^'(Pi) - ^-P-^^ < (7-15) 

C — 1 

^+{g{pi) - g{p+)) + ^+{p+ - pi)g\pi) + ^+p+^^ — < o . (7.i6) 

From the hypothesis (7.11) we can further reduce (7.15)-(7.16) to 

-(^(Pi) - g{p-)) + (Pi - P-)^'(Pi) > ^^P^ (7-17) 

{g{p+) - g{pi)) - (p+ - Pi)^'(pi) > ^^P+ • (7-18) 

Moreover, we observe from (7.1)-(7.4) that 

'^+(P+-Pi) = -v-{pi - P-)- 
Hence, in view of (7.11)-(7.12), we must have 

p_<pi<p+. (7.19) 

Let us note that 

{g{(^) - g{s)) -{(T- s)g\s) =11 g'\r)drdT 





r^'c-'** 


s 


/. r 


r 


/>'('■)*.. 


s 


/r I- 
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for every s < a. On the other hand, by simple algebra, we can compute g"{r) = p'{r)/r. 
Hence, the following equalities hold for every s < a: 

{9{a)-g{s))-{a-s)g'{s) 

and 

{g{s)-g{a)) + {a-s)g'{a)- 

As a consequence, and in view of (7.19), we can rewrite (7.17) and (7.18) equivalently as 

£'■ l" m,r,r > ^p_ , ,7.20) 

'-rm,r,r>^,,. ,7.21) 

Now, we introduce two new variables g_ and g+ defined by 

q- ■=p{pi) -p{p^), 

(1+ ■=P{P+) -P{pi)- 
Proving Lemma 7.2 is then equivalent to showing the existence of a pressure law p satisfying 
p{p+) — p{pi) = q+, p{pi) — p{p-) = 1- and for which the inequalities (7.20)-(7.21) hold. 
First, introducing / := p', we define the set of functions 

£:=|/gC°°(]0,oo[,]0,ooD : T / = ?- and ['^ f = q+ 
L J p- Jpl 

and the two functionals defined on C 

•Jpl J pi ^ 



L~{f):= p fy^drdr. 



' p- 

Therefore, a sufficient condition for finding a pressure function p with the above properties 
is that 

/+:=supL+(/)>^^p+ 

and 

/- := supL-(/) > ^- — p_. 
/e-c 2 

Let us generalize the space C as follows. We introduce 

A^^ := {positive Radon measures p on [pi,p+] : p([pi,p+]) = g+} , 
A^^ := {positive Radon measures p on [p_,pi] : p([p_,pi]) = g_} . 
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For consistency, we extend the functionals L^ and L~ defined on C to new functionals L. 
and L_ respectively defined on M.^ and on Ai~: 

L+(/i) := / / -dfi{r)dr for /x G M^ , 

•J pi J pi ^ 

rpi rpi I 

L^{fi) := / / -dfi{r)dT ioT fi e M^ ■ 



Upon setting 



and 



it is clear that 



771 "*" := max L+ia) 

IJ.&M+ 



m := max L_(yu), 



/"'' < rn^ and I < m 



Moreover, let us note that the the maxima m^m are achieved due to the compactness of 
Ai^ with respect to the weak* topology. By a simple Fubini type argument, we write 

L+{fi) = / — dp{r). 

J pi ^ 

Hence, defining the function h G C([pi, p+]) as h{r) := (p_|_ — r)/r allows us to express the 
action of the linear functional L+ as a duality pairing; more precisely we have: 

L+(yu) =< h,fi> for p G M~^. 

Analogously, if we define g G C([p_,pi]) as g{r) := (r — p^)/r, we can express L_ as a 
duality pairing as well: 

L^{fi) =< g,fi> for /i G Ai~. 

By standard functional analysis, we know that m^ must be achieved at the extreme points 
of Ai^. The extreme points of Ai^ are the single-point measures, i.e. weighted Dirac 
masses. For Ai^ the set of extreme points is then given by E^ := {q+S„ for a G [pi,p+]} 
while for Ai~ the set of extreme points is then given by E^ := {q-6„ for a G [p-,Pi]} In 
order to find m^, it is sufficient to find the maximum value of L± on E±. Clearly, we 
obtain 



and 





P+- Pl 
Pl 


', "-''] 


-q P'-P-. 



m = max 

o-e[pi,p+ 



TTi = max 

<^(^IP-,Pi] t' (T J ' Pl 



Furthermore, given the explicit form of the maximum points, it is rather easy to show that 
for every e > there exists a function f E C such that 

Pl 
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and 

r ff\ -^ P^~ P- 

L-U) > Q- ^• 

Pi 

Such a function / is the derivative of the desired pressure function p. D 

7.2. Part II of the proof of Lemma 7.1. We now choose pi = 1. Applying Lemma 7.2 

we set q± := ±{p{p±) — p{pi)) = ±{p{p±) — p(l)) and hence reduce our problem to finding 
real numbers p±, h'±, q±, a, (3, 7, 6, Ci satisfiyng 

i/_ < < Z/+ , < p_ < 1 < p+ , g± > (7.22) 

z/_(l-p_) = /3 (7.23) 

z/_(p_+«) =6 (7.24) 

Y-7 + g- = ^-/3 (7.25) 

z/+(l-p+) = /3 (7.26) 

u^{a-p^) = 6 (7.27) 

Y - 7 - 9+ = ^+/3 , (7.28) 

qM-p-) >^^P^ (7.29) 

g+(p+ - 1) > ^^P+ (7-30) 

and (5.11)-(5.12). 

Next, using (7.22), (7.23) and (7.26) we rewrite (7.29)-(7.30) as 

-Pq- > ^^ (-^-P-) (7.31) 

-/3g+ >^^z/+p+. (7.32) 

In order to simplify our computations we then introduce the new variables 

_ _ (^ 

/3 = —f3 , 6 = —6 , C = — , u^ = — z/_ , r+ = p+z/+ and r_ = p-i'^ = —p-U^ . (7.33) 

Therefore, our conditions become 

g±,r±,z/+,z/- >0 (7.34) 
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11)-(5.12) 


become 

a' + f 


< 


2C 











(7.35) 
(7.36) 
(7.37) 
(7.38) 
(7.39) 
(7.40) 

(7.41) 

(7.42) 



(7.43) 

(C-a^ + 7)(C - ^^ - 7) - (5 - aj)^ > . (7.44) 

We assume a^ 7^ 1 and solve for z/_, z/+ and r± in (7.35)-(7.38) to achieve 

_ 5 + ]3 S-]3 ~5-a(3 

V = -- — , Z/+ = and r± = . (7.45) 

1 + a 1 — a 1=fq; 

Observe that _ _ 

(5 - aPf 

r+r^ = — — . 

1 — a^ 

Hence, if we assume a^ < 1 and 6 > (3 > 0, we see that the z/"*", z/_,r-|- as defined in 
the formulas (7.45) satisfy the inequalities in (7.34). Hence, inserting (7.45) we look for 
solutions of the set of identities and inequalities 

a^ <1, 5>]3>0, q±>0 (7.46) 

C - 7 + g- = ^^ (7.47) 

1 + a 

C-^-q^ = J-fI--p (7.48) 

1 — a 

combined with (7.43) and (7.44). Observe that, if we assume in addition that C > ^, then 
a^ < 1, ~5 > (3 > Q and (7.49)-(7.50) imply the positivity of q±. We can therefore solve 
the equations (7.47)-(7.48) for q± and insert the corresponding values in the remaining 
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a^ 


< 1, (5>/3>0, 


?'>^ 
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1 + a 


> C- 

v 


'2J 


5 -a/3 
1 + a 


^ 


_ 1 - Q; 


> c- 
V 


1^ 

"2y 


5 -a/3 
1-a 


a^ 


+ ^^ < 2C 









inequalities (7.49) and (7.50). Summarizing, we are looking for a,/3,7,(5, C satisfying the 
following inequalities 



(7.51) 
(7.52) 

(7.53) 

(7.54) 
(7.55) 

(7.56) 

(7.57) 

(7.58) 

(7.59) 
(7.60) 

(7.61) 
(7.62) 
(7.63) 

(7.64) 

(7.65) 
(7.66) 



(C - a^ + 7)(C - /3 - 7) - (5 - al3f > . 
We next introduce the variable X = 6 — a/3 and rewrite our inequalities as 

a^ < 1, A> (l-a)/3>0, C> - 

^(l + a)(^'-C + 7)> (c-l^'-l)>^ 

l3{l-a)i-f + C-^)> fc-f-^)\ 

a^ + P^ < 2C 

{C-a^ + ^)(C-]3^-^)>\\ 

Observe that, if we require a, /3, 7 and C to satisfy the following inequalities 



a^ < 1 , C> 



1 



C - a^ + 7 > 
C - ^^ - 7 > 



(C - a2 + 7)(C - /3 - 7) > (1 - a)/3 > 
1 



P 



C 



C-a2 + 7>/3(l + a)VC-/3 -7 



then setting 



A:= V(C-a2 + 7)(C-/3 -^)-v, 
the inequalities (7.56)-(7.60) are satisfied whenever 77 is a sufficiently small positive number. 
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Observe next that (7.64) is surely satisfied if the remaining inequahties are and hence 
we can drop it. Moreover, if /3,7 and C satisfy 

^>0,C>^ (7.67) 

C-^^-7>0 (7.68) 

C-l+7>0 (7.69) 



/?' + ^ - CJ V^C-1 + 7 > 2/3 VC - /?' - 7 (7.70) 

then setting a = 1 — i), the inequahties (7.61)-(7.66) hold provided -i? > is chosen small 
enough. 

Finally, choosing C = |/3 , 7 = — 1/3 and imposing /3 > w | we see that (7.67), (7.68) 

and (7.69) are automatically satisfied. Whereas (7.70) is equivalent to 




5 ^/5' 

However the latter inequality is surely satisfied for /3 large enough. 

8. Classical solutions of the Riemann problem 

We show here that, if we restrict our attention to BV selfsimilar solutions of (1.1)-(1.3) 
which do not depend on the variable Xi, then for the initial data of Theorem 1.1 and 
Corollary 1.2 the solutions of the Cauchy problem are unique. We mostly exploit classical 
results about the 1-dimensional Riemann problem for hyperbolic system of conservation 
laws. We however complement them with some recent results in the theory of transport 
equations: the resulting argument is then shorter and moreover yields uniqueness under 
milder assumptions (see Remark 8.2 below). 

Proposition 8.1. Consider p{p) = p^ and any initial data of type (1.3) as in Lemma 
6.1 (which therefore include the data of the proof of Theorem 1.1 and Corollary 1.2). 
Then there exists a unique admissible self-similar hounded BVioc solution (i.e. of the form 
{p,v){x,t) = (r,u')(Y-)j of (1.1) with p bounded away from 0. 

Remark 8.2. In fact our proof of Proposition 8.1 has a stronger outcome. In particular 
the same uniqueness conclusion holds under the following more general assumptions: 

• p satisfies the usual "hyperbolicity assumption" p' > and the "genuinely nonlin- 
earity condition" 2p'{r) + rp"{r) > Vr > 0; 

• (p, v) is a bounded admissible solution with density bounded away from zero, 
whereas the BV regularity and the self-similarity hypotheses are assumed only 
for p and the second component of the velocity v. 

Remark 8.3. The arguments given below can be adapted to show the same uniqueness 
statement for the Cauchy problem corresponding to the data generated by Lemma 7.1. 
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This would only require some lengthier ad hoc analysis of the classical Riemann problem 
for the system (8.3), with p^ replaced by the pressures of Lemma 7.1. 

Proof. Observe that the initial data for the first component Vi is the constant — —. On 
the other hand: 

• p is a bounded function of locally bounded variation; 

• The vector field v = (0,f2) is bounded, has locally bounded variation and solves 
the continuity equation 

dtp + dw^{pv) = 0; (8.1) 

• f 1 is an L°° weak solution of the transport equation 

dtipvi) + div(p-i;t>i) = 

(8.2) 

-1(0,-) = -^. 

Therefore, the vector field v is nearly incompressible in the sense of [9, Definition 3.6]. 
By the BV regularity of p and v we can apply Ambrosio's renormalization theorem [9, 
Theorem 4.1] and hence use [9, Lemma 5.10] to infer from (8.1) that the pair (p, v) has the 
renormalization property of [9, Definition 3.9]. Thus we can apply [9, Corollary 3.14] to 
infer that there is a unique bounded weak solution of (8.2). Since the constant function is 
a solution, we therefore conclude that Vi is identically equal to — —. 

Set now m(x2,t) := p{x2,t)v2{x2,t). The pair p,m is then a self-similar BVioc weak 
solution of the 2x2 one-dimensional system of conservation laws 

dtp + dx^m = 

/ . X (8.3) 

dtm + d.,[^+p')=0, 

that is the standard system of isentropic Euler in Eulerian coordinates with a particular 
polytropic pressure. It is well-known that such system is genuinely nonlinear in the sense of 
[7, Definition 7.5.1] and therefore, following the discusssion of [7, Section 9.1] we conclude 
that the functions (p, m) result from "patching" rarefaction waves and shocks connecting 
constant states, i.e. they are classical solutions of the so-called Riemann problem in the 
sense of [7, Section 9.3]. It is well-known that in the special case of (8.3) the latter property 
and the admissibility condition determine uniquely the functions (p, rn). For instance, one 
can apply [17, Theorem 3.2]. D 
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